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EXECUTIVE  SUMMARY 


OBJECTIVE 

This  report  describes  an  approach  taken  to  test  the  practicality  of  an  instrument  that  uses  the 
reflection  coefficient  modulated  radiance  (shape  from  reflection)  to  measure  the  roughness  of  the 
ocean  surface.  Conceptually,  this  instrument  would  be  used  remotely  aboard  ship  to  sense  the  ocean 
wave  spectrum  of  the  ocean  surface;  ideally,  the  ocean  wave  spectrum  covered  by  this  instrument 
would  range  from  very  long  wavelengths,  on  the  order  of  decameters,  to  very  short  wavelengths,  on 
the  order  of  centimeters. 

RESULTS 

A  series  of  wave-tank  measurements  made  at  the  Scripps  Institution  of  Oceanography  (SIO), 
San  Diego,  CA,  clearly  show  that  the  technique  of  using  the  reflection  coefficient  modulated 
radiance  to  measure  surface  roughness  is  impractical.  The  technical  difficulties  are  associated  with 
the  failure  of  several  necessary  assumptions  of  shape  from  reflection.  The  first  assumption  is  that  the 
“sky”  illumination  source  is  nearly  isotropic;  in  practice,  this  assumption  is  not  valid.  The  second 
assumption  is  that  the  nonlinear  relationship  between  the  reflection  coefficient  and  wave  slope  is  not 
a  fundamental  limiting  factor;  in  practice,  the  sensitivity  of  the  reflection  coefficient  to  incidence 
angle  appears  to  be  a  limiting  factor. 

RECOMMENDATIONS 

It  is  recommended  that  the  technique  of  shape  from  reflection  be  abandoned  as  a  method  to  assess 
the  ocean  surface  roughness. 
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INTRODUCTION 


The  sea  surface  is  of  fundamentai  importance  to  U.S.  Navy  operations  not  onJy  because  it  is  the 
milieu  of  the  surface  fleet,  but  also  because  it  significantly  contributes  to  radar  clutter,  acoustic  and 
electromagnetic  signal  propagation,  and  boundary  layer  interactions.  Heavy  seas,  caused  by  a 
hurricane  or  a  typhoon,  may  extensively  modify  fleet  deployment  plans.  A  moderately  rough 
surface  may  increase  the  ambient  acoustic  noise  level,  thereby  reducing  the  maximum  detection 
range  of  underwater  targets.  Above  the  surface,  a  rough  sea  may  reduce  a  radar  system’s  detection 
range  by  decreasing  the  magnitude  of  the  reflection  coefficient. 

Statistical  measures  of  the  sea  surface  are  regarded  in  ocean  wave  analysis  as  the  only  descriptors 
that  are  significant,  both  observationally  and  theoretically,  because  the  surface  shape  varies  in  an 
irregular  and  unpredictable  manner  in  both  space  and  time.  These  statistical  measures  vary  in 
complexity  from  visual  estimates  of  the  significant  wave  height  to  measurements  of  the  energy 
spectra,  either  of  the  surface  function  (the  displacement  of  the  surface  from  some  reference  height), 
or  of  the  wave  slope.  The  distinctions  between  the  statistical  measures  are  related  to  the 
wavelengths-of-interest  and  to  the  instruments  used  in  the  measurements.  Long  wavelengths,  on  the 
order  of  decameters  and  greater,  affect  ship  maneuverability  and  are  associated  with  surface 
displacements  that  can  be  measured  by  satellite-borne  radars.  Moderate  wavelengths,  on  the  order  of 
meters  and  greater,  affect  acoustic  signal  propagation  and  can  be  measured  by  instruments  such  as 
buoy-mounted  accelerometers  and  resistive  wire  gauges  (wave  staffs).  Short  wavelengths,  on  the 
order  of  centimeters  and  less,  affect  electromagnetic  signal  propagation  and  are  associated  with 
small  surface  displacements  and  large  slopes  that  can  only  be  measured  using  noncontact 
electro-optic  techniques. 

High  sea  states  are  associated  with  long  wavelengths  and  are  generally  easy  to  predict  because 
these  conditions  are  associated  with  large-scale  forcing  functions  (hurricanes  or  typhoons)  that  are 
themselves  easy  to  detect  and  track.  The  time  and  length  scales  of  the  driving  forces  are  synoptic; 
that  is,  lengths  are  on  the  order  of  100  km  and  greater,  with  time  scales  on  the  order  of  days  to  weeks. 
Wind-generated,  short  wavelength,  small-gravity  and  capillary  waves  are  created  by  localized 
meteorologic  and  oceanic  conditions  that  are  difficult  to  predict  without  in  situ  measurements.  The 
scales  of  short  wavelength  driving  forces  are  on  the  order  of  a  kilometer  or  less,  with  time  scales  on 
the  order  of  seconds  to  minutes  (microscale).  Moderate  wavelength  seas  result  from  a  combination 
of  synoptic  scale,  mesoscale  (1  to  100  km,  minutes  to  days),  and  microscale  atmospheric  and  oceanic 
conditions. 

Onboard  Tactical  Decision  Aids  (TDAs)  have  been  developed  to  assist  operational  forces  in 
exploiting  or  mitigating  environmental  effects.  Ship  response  TDAs  can  be  used  to  assess 
maneuverability  in  various  sea  states.  Acoustic  and  electromagnetic  TDAs  assist  in  determining 
maximum  sensor  detection  ranges  against  selected  targets.  Although  computerized  tools  exist  to 
assess  the  effect  of  surface  roughness,  it  is  rarely  measured  directly,  but  is  usually  inferred  from 
observations  of  wind  speed,  or  from  estimates  of  the  significant  wave  height.  In  most  cases,  the 
infened  information  is  adequate  for  TDAs  because  they  were  developed  with  the  knowledge  that 
only  secondary  information  (e.g.,  wind  speed)  would  be  available  on  a  routine  basis.  If  additional 
information,  such  as  the  wave  spectrum,  were  to  become  readily  available,  most  of  the  TDAs  could 
be  updated  to  include  a  wave  spectrum  description  of  surface  roughness,  perhaps  substantially 
improving  the  prediction. 
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Although  a  number  of  instruments  exist  to  measure  the  spectra  of  wind-generated,  small-gravity, 
capillary  waves  in  laboratory,  or  near-laboratory  conditions,  currently,  none  are  available  for  use 
aboard  operational  ships.  Instruments  to  measure  the  moderate  wavelength  spectra  are  also  not 
applicable  for  use  aboard  operational  ships.  A  field  of  buoys  could  be  deployed  within  a  day  (at 
best),  but  is  generally  considered  impractical  because  of  the  expense,  and  because  the  buoy  field 
would  be  essentially  stationary.  Resistive. wire  gauges  are  too  delicate  for  shipboard  use,  and  they 
require  periodic  calibrations  and  cleaning. 

The  advantages  of  an  onboard,  in  situ,  wave  spectra  measurement  instrument  extend  beyond 
possible  improvements  in  TDA  prediction  capabilities.  The  improvements  in  basic  research  could 
be  considerable;  for  example,  such  an  instrument  would  valuable  to  both  acoustic  and 
electromagnetic  propagation  measurement  programs  to  firmly  establish  the  underlying 
relationships  of  surface  roughness  to  propagation.  Routine  surface  spectral  measurements  would 
benefit  both  global  climatology  and  pollution  studies. 

The  concept  of  such  an  instrument  is  straight  forward.  A  detector  (figure  1)  looks  down  on  the 
surface  of  the  water  and  measures  the  radiance  emitted  and  reflected  from  the  surface.  The  measured 
radiance  L  (in  two  dimensions  and  neglecting  any  path  transmittance  effects)  is 

L  =  (1/n)  1 1  [p.(l  -  0)  +  e.  p.  0]  dX  de  (1) 

where  Ps  is  the  black-body  radiant  emittance  of  the  surface;  q  is  the  Fresnel  reflection  coefficient 
evaluated  at  the  reflection  angle  Q;  Ea  is  the  emissivity  of  the  air  source;  Pa  is  the  black-body  radiant 
emittance  of  the  air  source;  X.  is  the  wavelength  of  the  radiation,  and  6  is  the  depression  angle  between 
the  detector  and  the  surface  element.  The  1/n  term  is  needed  for  Lambertian  surfaces.  If  Ps  and  Ea  Pa 
are  approximately  constant  over  X  and  6,  then  the  measured  radiance  L  is  modulated  by  q.  By 
selecting  X,  q  becomes  a  strong  function  of  Q,  which  is  itself  a  strong  function  of  the  surface  slope 
angle  a;  therefore,  the  radiance  can  be  strongly  connected  to  the  slope  of  the  wave.  If  the  spectrum 
of  the  wave  slope  is  known,  then  the  spectrum  of  the  surface  function  g  is  also  known. 


Figure  1.  The  geometry  of  an  instrument  to  measure  the  ocean  surface  height. 
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The  following  section  describes  an  approach  taken  to  test  the  practicality  of  an  instrument  that 
uses  the  reflection  coefficient  modulated  radiance  to  measure  surface  roughness.  Afta*  a  review  of 
the  physics  and  sensor  characteristics,  results  from  a  series  of  wave-tank  measurements  are 
examined.  It  is  concluded  that  the  technique  is  not  adequate  to  meastire  the  moderate-  and 
sh<xt- wavelength  spectrum.  The  nonlinear  nature  of  the  Q-a  relationship  and  the  assumption  of  a 
constant  sky  radiatix'  are  the  majw  sources  of  errm'.  Even  though  the  technique  is  not  satisfactrxy,  the 
methods  and  results  are  documented  not  only  to  prevent  future  studies  from  pursuing  this  path,  but  to 
encourage  and  challenge  other  investigator’s  to  examine  new  techniques  of  remotely  sensing  the 
ocean  surface  roughness  from  in  situ  measurements. 

SHAPE  FROM  REFLECTION 


The  technique  to  determine  the  surface  roughness  from  the  reflection  coefficient  modulated 
radiance  is  referred  to  as  ”sh^  from  reflection,”  (Jahne,  Waas,  and  Klinke,  1992).  Stilwell  (1969) 
was  the  Erst  to  prrqrose  using  photogrrqjhic  techniques  to  determine  the  surface  spectra  by  shiq)e 
from  reflection.  His  method  is  known  as  Stilwell  photography.  The  critical  assumptions  are  an 
isotropic  (ot  nearly  isotropic)  slqrlight  luminance  and  a  uniform  surface  temperature.  Then,  the 
radiance  depends  only  on  the  reflection  coefficient. 

Radiant  emittance  px  (mW  cm*^  (un*^)  is  related  to  black-body  temperature  and  the  wavelengtii 
of  the  radiation  through  the  Planck  equation; 

=  2  jc  cV[10X5(e*“/^'r  -  1)]  (2) 

where  c  is  the  speed  of  light;  h  is  Planck’s  constant:  Iq,  is  Boltzman’s  constant,  and  T  is  the 
temperature  (Kelvin). 

A  Lambotian  surface  is  a  perfectly  diffuse  surface  that  has  a  constant  radiance  independent  of 
viewing  direction.  The  constant  radiance  is  the  radiant  emittance  per  unit  solid  angle 

k  =  (3) 

which  is  expressed  in  units  of  mW  cm’^  sr^  pm'^. 

The  sensing  system  was  chosen  to  q)erate  in  the  far  infrared  for  Seva's!  reasons;  first,  in  the  8-  to 
12-micron  (pm)  band,  the  small  optical  depth  of  wato*  reduces  the  effects  of  internal  scatter 
(radiation  penetrating  the  surface  and  scattering  back  out  at  some  other  angle);  second,  unlike  some 
methods  thatrely  on  the  sun  as  a  source  of  photons,  the  sensors  would  be  able  to  operate  day  and  night 
(but  with  fewer  available  photons);  third,  sensors  operating  in  the  8-  to  12-pm  band  were  available 
for  use.  These  sensors  were  Everest  Interscience  Inc.,  Model  4000,  that  had  been  purchased  to 
support  other  experimental  efforts.  The  ffeld-of-view  (fov)  angle  for  these  sensors  is  4  degrees. 

Figure  2  shows  the  black-body  temperature  with  respect  to  the  band  averaged  8-  to  12-pm 
radiance  L.  The  T-L  relation  is  slightly  nonlinear  over  the  entire  temperature  range  shown,  from  0  to 
40  degrees;  however,  over  smaller  temperature  segments  (of  about  10  degrees),  the  T-L  relation  is 
nearly  linear. 

In  the  infrared,  the  refractive  index  of  water  q  is  well  known  (Hale  and  (Juerry,  1973).  Figure  3 
shows  the  relationship  between  the  complex  refractive  index  of  saline  water  and  wavelength  (values 
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Figure  2.  Band  average  black-body  radiance  in  relation  to  the 
black-body  tempotiture. 
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Figure  3.  The  refractive  index  of  saline  water.  The  values 
correspond  to  water  taken  from  the  Pacific  Ocean  (Querry  et  al., 
1977). 
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fr(Hn  Querry  et  ol.,  1977,  table  3,  Pacific  Ocean);  the  circles  rq>resent  the  real  part  of  t),  and  the 
diamonds  rqniesent  die  imaginary  part  of  t].  The  differences  between  fresh  water  and  saline  water  t] 
are  insignificant  for  the  current  analysis. 

At  an  air-water  boundary,  for  3  specific  radiation  wavelength,  the  Fresnel  reflection  coefficient, 
QX  is  a  fimctitHi  of  die  reflecdo’'  angle,  fi,  and  the  cmnplex  refractive  index  of  water,  t),  accoiding  to 

1/2 

a  =  |4  m?m?  +  [m,  -  mf  -  sin^(Q)]^| 
p^  =  1/2  [-  m?  +  m?  +  sin^(Q)  +  a] 

=  1/2  [m?  -  m?  +  sin^(Q)  +  a] 

(q  -  cos(Q))2  +  p^ 

- 

(q  +  cos(Q))^  +  p^ 

[(m,  -  m?)  cos(Q)  —  q]^  +  (2  mr  m^  cos(Q)  —  p)^ 

[(m,  —  m?)  cos(Q)  +  q]^  +  (2  m,  mj  cos(Q)  -I-  p)^ 

ex=  1/2[Qvx  +  Qwl]  (4) 

where  mr  is  the  real  part  of  q;  mi  is  the  imaginary  part  of  q;  Q),x  is  the  reflection  coefficient  for 
polarization  perpendicular  to  the  plane  of  incidence  (h(»izontal  polarization);  QvX  is  die  reflection 
coefficient  fcx  polarization  parallel  to  die  plane  of  incidence  (vertical  polarization);  qx  is  the 
reflection  coefficient  fcx  impolarized  light;  and  the  refractive  index  of  air  is  assumed  to  be  one 
(Stratton,  1941).  Figure  4  shows  the  magnitude  of  the  band-averaged,  unpolarized,  Fresnel 
reflection  coefficient  q  with  relation  to  the  reflection  angle  Q.  Fm*  angles  near  grazing,  the  surface  is 
a  reflects;  whereas,  fcx*  angles  near  ntxmal  incidence,  die  surface  is  an  emitter. 

The  nonlinear  f<xtn  of  the  band-avo-aged  q,  in  terms  of  the  wave-slope  angle  a,  is  shown  in 
figure  S  for  a  sensw  depression  angle  0  of  20  degrees.  The  magnitude  of  q  changes  slightly  fm 
positive  wave  slc^ies  and  changes  dramatically  for  negative  wave  slqies.  This  implies  that  (1)  the 
sensOT  should  not  measure  a  significant  change  in  radiance  when  it  is  looking  at  the  front  of  a  wave 
(assuming  that  the  wave  is  traveling  toward  the  senstH*),  and  (2)  the  sensor  should  measure  a  large 
change  in  radiance  when  it  is  looking  at  the  back  side  of  a  wave.  The  nonlinear  Q-a  relationship 
should  cause  troubles  in  determining  the  surface  spectra,  particularly  if  the  sensor  has  a  narrow 
field-of-view  angle.  Additional  complications  arise  fin*  cases  of  negative-angle  slopes  (back  side  of 
a  wave),  where  there  is  the  possibility  of  shadowing  and  double  surface  reflections. 

Although  thecxy  clearly  indicates  considerable  difficulties  with  using  shape  from  reflection,  a 
set  of  field  measurements  provided  enough  tantalizing  infixination  to  continue  pursuing  the  effort. 
In  1991,  the  Everest  Interscience  sensors  were  being  used  to  monitor  the  sea-sinf^ace  temperature  to 
support  a  study  of  radar  detection  of  low-altitude  targets.  These  IR  sensors  were  positioned  22  m 
above  the  ocean  surface.  Unfrxtunately,  the  senscx  fov  was  restricted  to  looking  in  the  surf  zone, 
where  breaking  waves,  sea  foam,  and  other  contaminants  were  present;  however,  with  a  slight 
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Figure  4.  The  band  average  Fresnel  reflection  coefficient  for 
unpolarized  light  at  an  air-sea  interface  in  relation  to  the 
reflection  angle  Q. 
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Figure  5.  The  band  average  Fresnel  reflection  coefficient  for 
unpolarized  light  at  an  air-sea  interface  in  relation  to  the  wave 
slope  angle  a.  The  sensor  depression  angle  6  is  20  degrees. 
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modification  to  the  existing  sensor  system,  a  usable  IR  picture  of  the  surf  zone  was  made.  Figure  6 
shows  a  trend-removed  time  history  of  the  IR  signature  in  terms  of  temperature.  The  wavelike  nature 
of  the  surf  zone  is  clearly  represented  by  the  temperature  fluctuations.  First  impressions  of  the  data 
seem  to  indicate  that  changes  in  temperature  amplitude  could  be  related  to  the  wave  height;  for 
example,  figure  6  seems  to  indicate  a  large  wave  followed  by  three  (or  four)  smaller  waves,  followed 
by  three  moderate  waves,  followed  by  a  small  wave,  which  is  followed  by  three  larger  waves.  This 
pattern  of  large  and  small  waves  could  not  be  easily  verified  by  a  separate  measurement.  At  the  time, 
however,  the  wave  structure  in  the  surf  zone  was  observed  visually,  and  the  height  observations 
generally  agree  with  the  measured  data;  that  is,  larger  waves  appeared  to  yield  larger  temperature 
changes.  A  visual  observation  of  the  swell  indicated  a  wave  period  of  about  15  seconds.  The 
periodicity  is  evident  from  the  measured  temperature  data. 
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Figure  6.  An  IR  signature  of  waves  in  the  surf  zone. 


It  is  not  possible  to  link  the  surf  zone  measurements  directly  to  theory  because  of  the  surface 
contaminants  (foam,  sea  weed,  etc.)  and  their  unknown  effect  on  the  radiance  (temperature); 
however,  the  results  were  enticing  enough  to  warrant  testing  the  system  under  laboratory  conditions. 
Funding  was  made  available  in  January  1992,  and  a  series  of  wave-tank  measurements  were  made. 
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WAVE-TANK  MEASUREMENTS 


EQUIPMENT  CONFIGURATION 

Figure  7  illustrates  Ae  equipment  set  up  at  the  \^d-Wave  Channel  facility  at  the  Scripps 
Institution  of  Oceanography  (SIO),  San  Diego,  California.  Interim-  dimensions  of  the  tank  are 
44.S-m  long,  2.39-m  wi^,  and  2.44*m  deep  with  a  maximum  water  depth  of  2  m.  The  wave 
generatm  is  an  electro-hydraulic,  servo-controlled  paddle  that  can  create  a  maximum  wave  height  of 
60  cm.  A  U4)ered  beach  at  the  far  end  of  the  tank  attenuates  reflections,  and  the  glass- walled  tank  can 
be  enclosed  on  the  top  to  allow  blower-generated  wind  velocities  up  to  16  m/s  over  the  water  surface. 


Figure  7.  The  equipment  as  set  up  at  the  SIO  Wind- Wave 
Channel  facility  in  San  Diego,  CA.  The  IR  sensor  was 
1 12  cm  above  the  still,  fresh  water  surface  and  tilted 
downward  at  an  angle  of  20  degrees.  A  wave  staff  was 
positioned  at  the  still-surface,  beam-center,  specular 
reflection  point. 

Wind  effects  on  the  surface  were  not  examined.  Only  paddle  generated  waves  wme  measured. 
The  tank  was  Hlled  to  a  depth  of  145  cm  with  fresh  water.  Diere  may  have  been  some  contamination 
of  the  surface  because  it  was  not  skimmed  before,  or  during,  the  3  days  of  measurements  (2-4  June, 
1992);  however,  it  was  felt  that  if  the  technique  did  not  work  with  a  fairly  clean  facility  there  would 
be  no  chance  of  it  working  in  the  freld. 

SIO  provided  the  wave  staffs  with  their  associated  oscUlatm  and  buffer  amplifier  circuits.  The 
signal  ouq)ut  was  calibrated  for  wave  height  by  using  a  graduated  scale  attached  to  the  wave  staff. 
While  the  staff  was  moved  up  and  down  in  sdll  water,  the  height  above  the  surface  and  the  signal 
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voltage  were  recorded.  Figure  8  shows  a  calibration  curve  that  was  recorded  on  3  June  1 992  at  1 230 
PDT,  the  circles  represent  the  measured  points.  The  ordinate  is  the  surface  function  or  surface  height 

^  (in  cm),  which  is  measured  relative  to  the  still  surface.  A  linear  fit  to  this  calibration  curve  yields  a 
cmrelation  coefficient  of  0.999969.  Most  of  die  other  calibrations  made  during  the  measurement 
period  were  also  nearly  linear.  Although  not  direcdy  measured,  schematics  of  the  wave  staff  elec¬ 
tronics  indicate  that  the  high  frequency  cutoff  of  the  output  amplifier  circuit  is  about  6  Hz.  The  rol¬ 
loff  of  the  response  is  estimated  to  be  about  6  dB  per  octave. 


SENSOR  VOLTAGE 


Figure  8.  A  typical  wave  staff  voltage-versus-surface  height 
calibration. 

PriOT  to  installation  at  the  SIO  facility,  an  Everest  Interscience,  Model  4(X)0,  IR  sensor  and  its 
signal  conditioning  circuitry  were  tested  to  determine  the  frequency  response.  A  camera  shuttu-, 
with  an  qierture  of  1.4  cm,  was  placed  between  the  IR  sensm:  and  a  heated  pot  of  water.  The  condi¬ 
tioned  output  voltage  was  recorded  prior  to,  during,  and  after  opening  the  camera  shutter.  Figure  9 
shows  a  typical  response  curve.  The  rise  time  of  the  leading  edge  is  about  12.5  ms  (approximately  2 
ms  is  due  to  the  shutter);  the  effective  time  constant  of  the  system  is  about  33  ms. 

It  was  decided  to  set  the  analysis  cutoff  frequency  at  10  Hz,  a  fiequency  greater  than  the  high 
frequency  cutoff  of  the  wave  staffs  and  about  one-third  of  the  IR  transducer  high  frequency  cutoff. 
IVo  ind^ndent  11-pole  Chebyschev,  low-pass  filters  were  constructed  using  Datel  (Mo^l  FU- 
D5LA1  and  FLJ-D6LA1)  active  filters.  The  voltage- versus-frequency  response  of  tiie  anti-aliasing 
filters  is  shown  in  figure  10.  At  a  frequency  of  10  Hz,  the  voltage  response  is  2  dB  down;  at  15  Hz, 
the  voltage  response  is  greater  than  30  dB  down. 
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TEMPERATURE  (°C) 


An  IR  senscH*  head  was  positioned  in  the  center  of  die  Mnd-Wave  Channel  at  1 12  cm  above  the 
still  surface  and  approxima^y  10  m  from  the  wave  generating  paddle.  The  dqiression  angle  0  of  the 
sensor  was  set  to  20  degrees  by  using  an  inclinometer.  A  resistive  wire  wave  staff  was  posititmed 
about  307.7  on  upstream  near  the  edge  of  the  channel  at  the  still-surface  beam-center  specular 
reflection  point  of  the  sensor  (figure  7).  Later  in  the  measurement  period,  a  second  wave  staff  was 
added  69.8  cm  downstream  from  the  first  to  determine  the  phase  velocity  of  the  waves. 

Analog  data  from  both  the  IR  and  the  wave  staff  were  digitized  by  a  12-bit  analpg-to-digital 
(A/D)  converter  that  sampled  the  data  at  a  32-Hz  rate.  The  A/D  converter,  a  National  Instruments, 
Model  MIO-16,  required  a  minimum  of  25-ps  delay  between  samples.  Compared  to  the  cycle 
interval  of  31.25  ms,  this  delay  is  negligible  and  is  igntxed  in  the  following  analysis.  All  data  were 
receded  under  control  of  an  IBM-compatible  microcomputer.  Typically,  1 1  sets  of  1024  wcx’ds  per 
channel  were  stored  per  recording  run. 

Figure  11  shows  the  deep-water  dispersion  relation  for  free  waves  on  still  water.  The  angular 
frequency  is 

0)2  =  g  k  -f-  72  (5) 

where  g  is  the  acceleration  of  gravity  (980  cm/s^),  and  k  is  tee  angular  wave  numbo*  Infk  (Lamb, 
1945).  F<n-  gravity  waves,  wavelengths  »  1.7  an,  tee  phase  velocity  is  found  to  be 

C  =  g/o).  (6) 

Figure  12  shows  tee  cross-ctxrelation  between  two  wave  staffs  that  were  spatially  sq)arated  by 
69.8  cm.  The  wave  staff  recwds  were  made  with  a  driving  frequency  of  1  Hz.  The  time  to  the  first 
positive  peak  is  0.4687  s  and  the  measured  phase  velocity  Co  is  148.9  an/s.  A  comparison  to  the 
calculate  phase  velocity  C  (156  cm/s)  indicates  good  agreement  between  deep-water  theory  and 
measurements.  Measured  phase  velocity  at  wave  frequencies  of  1 .0, 0.8,  and  0.6  Hz  are  shown  as  tee 
three  circles  in  figure  11. 

The  tests  were  conducted  at  tee  SIO  facility  over  a  three-day  period,  from  2-4  June  1992.  Much 
of  the  first  day  was  involved  in  equipment  set  up  and  preliminary  testing.  On  tee  second  day,  testing 
was  directed  toward  specular  reflection  measurements  of  a  strong  IR  source  (a  pcntable  electric 
heater)  and  ambient  measurements  of  different  driving  frequencies  with  tee  tank  covered  in  black 
plastic.  The  tank  was  covo-ed  to  reduce  stray  reflections  from  tee  building  interior.  Testing  on  the 
last  day  was  directed  toward  ambient  measurements  of  different  driving  frequencies  with  tee  tank 
open  to  the  ‘*sky”.  The  data  from  tee  4  Jime  measurement  poiod  is  examined  in  tee  following 
section. 

4  JUNE  1992  MEASUREMENTS 

On  4  June  1992,  a  series  of  wave-tank  measurements  wae  made  for  wave  frequencies  of  1 .0, 0.8, 
0.6,  and  0.4  Hz.  Figure  13  shows  tee  surface  height  g  and  wave  slope  angle  a  for  a  1-Hz  wave  of 
about  7-cm  amplitude  measured  at  1205  PDT  (1024  samples  taken  over  a  32-second  time  interval). 
The  abscissa,  wave  position,  is  the  spatial  coordinate  (in  cm)  teat  would  exist  if  tee  time  coordinate 
were  multiplied  by  the  phase  velocity  Cq.  A  wave  position  of  0  is  defined  as  tee  first  sample  in  tee 
data  reco'd.  Multiple  waves  are  obsoved;  wave  slopes  vary  from  a  maximum  of  about  +10  to  a 
minimum  of  about  - 10  degrees.  Figure  14  expands  the  20-  to  30-m  portion  of  the  wave  position  for 
greater  detail.  Wave  height  and  slope  are  shown  to  be  reasonably  well  behaved. 
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WAVELfNGTH  (cm) 

Figure  11.  The  dispersion  relation  for  free  waves  on  still  water.  The 
circles  indicate  measured  phase  velocities  at  frequencies  of  1.0, 0.8, 
and  0.6  Hz. 
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Figure  12.  The  cross  correlation  of  two  wave  staffs  separated  by 
69.8  cm.  The  measured  phase  velocity  for  a  driving  frequency 
of  1.0  Hz  is  148.91  cm/s. 
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WAVE  POSITION  (cm) 


Figure  13.  The  surface  height  ^  and  wave  slope  angle  a  for  a 
1-Hz  wave  of  iqjproximately  7-cni  amplitude. 


WAVE  POSITION  (cm) 

Figure  14.  A  detailed  view  of  the  20-  to  30-m  ptntion  of  the 
wave  position  shown  in  the  previous  figure.  Surface  height  and 
wave  slope  behave  as  expected. 
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Figure  15  expands  the  C  data  in  even  greater  detail  to  illustrate  the  method  of  shape  from  refrac¬ 
tion.  Each  circle  on  this  figure  represents  a  sample  of  the  surface  height  as  record^  by  the  wave 
staff.  Assume  that  at  some  instant  in  time,  say  tj,  the  wave  staff  measures  the  surface  height  at  the 
wave  position  indicated  by  the  arrow;  surface  heights  to  the  left  of  the  arrow  were  recorded  {xior  to 
time  ti;  whereas,  surface  heights  to  the  right  of  the  arrow  were  recorded  after  time  t;.  At  exactly  time 
ti,  die  IR  sensor’s  fov  (depicted  as  the  shaded  region  in  figure  15)  falls  on  1 1  surface  facets:  5  facets  to 
the  left  of  the  arrow  position,  and  6  facets  to  the  right  of  the  arrow  position.  For  each  of  these  11 
facets,  the  wave  slope  an^e  a  is  approximated  as 

a  =  tan"^(A^/Ax)  (V 

where  AC  is  the  height  difference  between  consecutive  samples,  and  Ax  is  the  position  difference 
given  by  Co/32.  The  depression  angle  6  is  calculated  by  knowledge  of  the  geometry  between  the 
center  of  the  facet  and  the  IR  sensor  (figure  7).  The  reflection  angle  Q  is  related  to  6  and  a  as 

Q  =  Jt/2  -  6  -  a  (8) 


WAVE  POSmON  (cm) 

Figure  15.  The  24-  to  25-m  portion  of  the  wave  position  shown  in  figure  13.  The 
shaded  area  represents  the  IR’s  sensor’s  field-of-view  (figure  7)  at  the  moment  ii 
time  when  the  wave  staff  measured  the  wave  height  indicated  by  the  arrow.  The 
reflectance  and  emittance  of  the  wave  facets  within  the  IR  sensor’s  fov  are 
summed  to  determine  the  temperature  of  the  surface  (shown  by  the  chain  line 
linking  the  diamond  symbols). 


The  Fresnel  reflection  coefficient  is  readily  found  using  the  Q-q  relation  shown  in  figure  4. 
'nut)ugh  the  T'L  connection  of  figure  2.  the  bulk  water  temperature  7$  and  the  “sky”  ambient  air 
temperature  Ta  that  are  recorded  during  the  measurement  period  determine  the  band-averaged  black- 
body  radiance  from  both  the  water  surface  and  tiie  air  source.  The  values  of  7$  and  Ta  for  this  mea¬ 
surement  are  20.9  and  2 1 .6  Celsius,  respectively.  The  band-averaged  ambient  emissivity  is  assumed 
to  be  one;  therefore,  the  radiance  as  measured  by  the  IR  sensm*  at  time  ti  is  approximated  as 

I  -  ej)  +  f.ej]  (9) 

j»i-n 

where  the  indices  m  and  n  are  chosen  to  satisfy  the  conditions  of 

112  -(307.7  -I- m  Co/32)  tan(18*)>  Si +  m,  and 

112  -  (307.7 -I- n  Co/32)  tan(22«)  <5i+n.  (10) 

These  conditions  simply  require  that  the  center  of  tiie  wave  facet  be  in  the  IR  sensw  field-of- 
view. 

The  radiance  Li  is  converted  to  a  tenq)erature  Ti  tiirough  the  T-L  connection  of  figure  2.  The 
surface  tmnperature  fluctuation  t  (in  units  of  centittegree  Celsius)  is  defined  as 

Ti  =  100  [Ti  -  (ao  +  ai  Xi)!  (11) 

where  Xi  is  the  wave  position  at  time  ti,  ao  and  ai  are  the  coefficients  of  a  linear  least  squares  fit  to  the 
T-x  pairs  calculated  for  the  data  recwd  (typically  1024  pairs). 

The  surface  temperature  fluctuation  calculated  fi'om  knowledge  of  ^ ,  T,,  Ta,  and  Co  is  shown  in 
figure  15  as  the  chain  curve  connecting  the  diamond  symbols.  At  time  q,  the  surface  tempoatuie 
fluctuation  is  about +1.3  centidegree  Celsius.  Ilie  shiqje  of  the  predicted  t  curve  is  a  distorted  sinu¬ 
soid;  its  negative  peak  occurs  slightly  before  the  negative  peak  of  C,  but  its  positive  peak  occurs 
considerably  befOTe  the  positive  peak  of  ^ .  The  raggedness  of  the  positive  peak  is  attributed  to  the 
distmtions  in  the  Q-a  relationship  as  a  transitions  from  positive  to  negative  angles  (figure  S).  Over 
the  20-  to  30-m  portion  of  the  wave  position  shown  in  figure  16,  the  predicted  x  curve  illustrates  the 
same  characteristics  as  the  detailed  view. 

Figure  17  compares  the  predicted  to  the  measured  surface  temperature  fluctuation.  The  mea¬ 
sured  curve  (solid)  appears  to  be  offset  in  phase  and  also  appears  to  be  modulated;  it  has  the  charac¬ 
teristics  of  a  beat  frequency.  The  modulation  is  evident  in  Ae  comparison  of  predicted  to  measured  x 
shown  in  figure  18.  Although  the  magnitude  of  the  predicted  and  measured  x  agree  reasonably  well 
(which  indicates  that  the  assumptions  used  in  shiq>e  from  reflection  themy  are  reasonable),  the  cor¬ 
relation  of  prediction  and  measiuement  are  disappointing. 

Figures  19-22  show  the  cross  correlation  between  the  measured  and  predicted  x  fw  driving  fre¬ 
quencies  of  1.0, 0.8, 0.6,  and  0.4  Hz,  respectively.  Each  lag  is  1/32  of  a  second.  Table  1  lists  the 
values  of  the  correlation  coefficient  and  the  lag  for  the  first  three  frequencies.  The  0.4-Hz  driven 
frequency  experienced  severe  distortions  in  C  that  are  probably  related  to  powly  attenuated  tank 
reflections  and  its  cmrelation  is  low  (figure  22).  The  correlation  fw  the  first  three  frequencies  (there 
is  ctmsiderable  confidence  in  these  measurements)  show  a  marginal  value  of  about  0.5.  An  indica¬ 
tion  of  problems  with  the  sh^  from  reflection  method  is  seen  by  the  sign  reversal  of  the 
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Figure  16.  The  surface  height  and  the  predicted  surface  tem> 
perature  fluctuation  (in  centidegree  Celsius)  for  the  20-  to 
30-ni  portion  of  the  wave  position  shown  in  figure  13.  The  pre 
dieted  surface  temperature  appears  as  a  distorted  sinusoid. 


WAVE  POSITION  (cm) 


Figure  17.  A  comparison  of  the  measured  (solid  line)  and 
predicted  surface  temperature  fluctuation.  The  measured 
temperature  appears  to  be  modulated;  whereas,  the  predicted 
temperature  is  not  modulated. 
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Figure  18.  An  expanded  view  of  the  predicted  and  measured 
surface  temperature  fluctuation.  The  beat  frequency  nature 
of  the  measured  temperature  is  readily  seen. 


UG  (MEASURED  t) 

Figure  19.  The  cross  correlation  of  the  measured  and 
predicted  temperature  fluctuation.  Driving  frequency  is 
1.0  Hz  at  an  amplitude  of  about  7  cm. 
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Figure  20.  The  cross  correlation  of  the  measured  and 
predicted  temperature  fluctuation.  Driving  frequency  is 
0.8  Hz  at  an  amplitude  of  about  10  cm. 
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Figure  21.  The  cross  correlation  of  the  measured  and 
predicted  temperature  fluctuation.  Driving  frequency  is 
0.6  Hz  at  an  amplitude  of  about  10  cm. 
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Figure  22.  The  cross  corelation  of  the  measured  and  pre¬ 
dicted  temperature  fluctuadtm.  Driving  frequency  is  0.4  Hz 
at  an  amplitude  of  about  8  cm. 

correlation  coefficient  between  driven  frequencies  of  1  and  0.6  Hz.  From  table  1,  the  first  peak  in 
correlation  is  about  0.5  for  die  1-  and  0.6-Hz  frequencies;  whereas,  for  the  0.8  Hz-frequency,  the  first 
wave  peak  in  correlation  is  about  -0.64. 


Table  1.  Correlation  between  measured  and  predicted  tenqierature  fluctuations. 
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Figure  23  shows  surface  height  C  and  wave  slope  a  curves,  and  figure  24  provides  measured  and 
predicted  t  curves  for  the  first  16  seconds  of  data  recorded  with  a  driving  firequency  of  0.4  Hz.  The 
nonsinusoidal  sluq)e  of  C  is  evident  (cf  figures  14  and  23).  The  measured  phase  velocity  Co  of  297.8 
cm/s  does  not  compare  to  the  predicted  phase  velocity  C  of  389.9  cm/s;  a  measurement  error  is  sus¬ 
pected;  however,  the  measured  phase  velocity  was  used  to  calculate  the  wave  sl(^)es .  The  slopes  are 

strongly  affected  by  the  C  distcations  and  have  considerable  impact  on  the  calculation  of  q  and  the 
predict^  t.  Even  ^ough  the  correlation  of  measured  and  predicted  t  is  very  pota*  (figure  22),  the 
compaiisons  shown  in  figure  24  are  qualitatively  good.  There  are  deviations  between  the  two 
curves,  but  thov  are  also  periods  of  agreement.  The  magnitudes  are  in  general  agreement,  and  the 
trends  appear  visually  related.  In  particular,  the  comparison  in  the  interval  from  10  to  14  seconds  is 
remarkable.  One  can  argue  die  fact  that  even  randmn  data  align  sometimes,  but  the  overall  sense  of 
agreement  of  the  measured  to  predicted  t  was  encouraging;  however,  the  intent  of  shape  from 
reflection  is  to  relate  the  measured  t  to  C .  The  results  of  these  comparisons  are  dis^pointing. 

Figures  25-28  are  the  cross  correlation  between  the  measured  t  and  the  measured  ^  for  wave 
frequencies  of  1 .0, 0.8, 0.6,  and  0.4  Hz  respectively.  Table  2  lists  the  values  of  the  correlation  coeffi¬ 
cient  and  the  lag  for  the  four  frequencies.  Even  though  the  0.4-Hz  wave  frequency  experienced 
severe  distortions  in  C  (figure  23),  the  aoss  ctxrelation  coefficient  is  not  substantially  different  from 
the  coefficient  for  the  other  frequencies;  however,  the  magnitude  of  the  ‘'average”  cmrelation  coeffi¬ 
cient,  about  0.6,  is  maiginal.  In  addition,  the  sign  change  of  the  ctxrelation  coefficient  between  fre¬ 
quencies  is  disturbing. 


Table  2.  Correlation  between  measured  temperature  fluctuations  and  surface  height. 
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Further  evidence  of  the  limitations  of  shape  from  reflection  is  shown  in  figures  29-32,  which  are 
the  distributions  of  measured  t  and  measured  C  fcx*  wave  frequencies  of  1.0, 0.8, 0.6,  and  0.4  Hz, 
respectively.  These  figures  are  derived  from  10,240 samples  (lOdatarecwdings  of  1024  samples)  of 
each  frequency.  The  height  distribution  is  about  what  one  would  e]q)ect  far  a  sinusoidal  function.  Its 
classic  cusp  shi^  is  expected  to  peak  near  maximum  and  minimum  values  of  C .  The  distribution  of 
measured  surface  temperature  t  has  no  relationship  to  the  ^  distribution.  If  anything,  thex  distribu- 
ticms  shown  in  figures  29-32  resemble  a  distribution  of  a  random  process,  not  a  sinusoidal  process. 
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Figure  23.  Surface  height  C  and  wave-slope  angle  a  for 
waves  of  8  cm  amplitude  at  a  driving  frequency  of  0.4  Hz. 
The  height  and  slope  are  not  sinusoidal. 
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Figure  24.  Measured  and  predicted  surface  temperature  fluc¬ 
tuation  for  waves  of  8  cm  amplitude  at  a  driving  frequency 
of  0.4  Hz. 
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Figure  25.  The  cross  correlation  of  the  measured  temperature 
fluctuation  and  the  measured  surface  height.  Driving  fre¬ 
quency  is  1.0  Hz  at  an  amplitude  of  about  7  cm. 
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Figure  26.  The  cross  correlation  of  the  measured  temperature 
fluctuation  and  the  measured  surface  height.  Driving  fre¬ 
quency  is  0.8  Hz  at  an  amplitude  of  about  10  cm. 
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Figure  27.  The  cross  correlation  of  the  measured  temperature 
fluctuation  and  the  measured  surface  height.  Driving  fre¬ 
quency  is  0.6  Hz  at  an  amplitude  of  about  10  cm. 


Figure  28.  The  cross  correlation  of  the  measured  temperature 
fluctuation  and  the  measured  surface  height.  Driving  fre¬ 
quency  is  0.4  Hz  at  an  amplitude  of  about  8  cm. 
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Figure  29.  The  distribudons  of  measured  surface  height  t 
and  measured  temperature  fluctuadon  x.  Driving  frequency 
is  1.0  Hz  at  an  amplitude  of  about  7  cm. 
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Figure  30.  The  distributions  of  measured  surface  height  ^ 
and  measured  temperature  fluctuation  x.  Driving  frequency 
is  0.8  Hz  at  an  amplitude  of  about  10  cm. 
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Figure  31.  The  distributions  of  measured  surface  height  t 
and  measured  temperature  fluctuation  t.  Driving  frequency 
is  0.6  Hz  at  an  amplitude  of  about  10  cm. 
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Figure  32.  The  distributions  of  measured  surface  height  C 
and  measured  temperature  fluctuation  x.  Driving  frequency 
is  0.4  Hz  at  an  amplitude  of  about  8  cm. 
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On  4  June  1992,  measurements  were  made  with  the  tank  open  to  the  sky.  On  the  previous  day, 
measurements  were  made  with  the  tank  covered  in  black  plastic  to  minimize  angular  variations  in  the 
sky  luminance  function.  However,  irrespective  of  whether  or  not  the  tank  was  covered,  there  is  no 
significant  difference  in  the  comparisons  of  measured  to  predicted  data;  the  4  June  measurements  are 
typical  of  what  was  seen  with  the  tank  covered.  In  the  following  paragraphs,  measurements  from 
both  days  are  used  to  examine  the  surface  height  variance. 

SURFACE  VARIANCE 

The  variance  of  the  surface  height  is  a  useful  and  readily  observable  statistic.  In  its  simplest 
form,  the  variance  is 

v2  =  i/(N  -  1)2]  [x  -  <  X  >]2  (12) 

where  N  is  the  number  of  x  variables,  and  <x>  is  the  average  or  expected  value  of  x.  In  a  more  com¬ 
plicated  form,  the  power  spectrum  is  the  frequency  distribution  of  variance. 

Figure  33  compares  the  variance  of  the  measured  surface  temperature  t  to  the  variance  of  the 
measured  surface  function  ^  for  four  driving  frequencies.  The  solid  (filled)  symbols  correspond  to 
data  taken  on  3  June  1992  (tank  covered);  the  open  (not  filled)  symbols  correspond  to  data  taken  on  4 
June  1992.  Each  symbol  (solid  or  open)  is  the  variance  of  1024  samples  (N  =  1024).  It  is  obvious  that 
the  surface  temperature  variance  is  not  consistently  related  to  the  variance  of  the  height. 
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Figure  33.  Comparison  of  measured  surface  temperature 
variance  and  measured  surface  height  variance.  The  open 
symbols  represent  data  that  were  taken  with  the  tank  open  to 
the  sky;  the  filled  symbols  represent  data  that  were  taken 
with  the  tank  covered. 
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CONCLUSIONS 


Shape  from  reflection  is  not  a  viable  technique  to  determine  the  simplest  statistics  required  to 
describe  the  ocean  surface.  Although  the  technique  appears  to  have  merit  (e.g.,  figures  6  and  24),  the 
qualitative  is  not  supported  by  the  quantitative.  Correlations  of  measured  and  predicted  temperature 
fluctuations  (figures  19-22  and  table  1)  are  poor;  the  average  correlation  coefficient  is  about  0.5. 
Correlations  of  measured  temperature  fluctuations  and  surface  height  (figures  25-28  and  table  2)  are 

*  also  poor,  averaging  about  0.6.  These  correlations  are  for  near  laboratory  conditions;  it  is  expected 
that  the  correlations  in  actual  field  conditions  would  not  be  any  better. 

The  death  knell  of  shape  from  reflection  is  provided  by  figures  29-32.  Hie  comparisons  of  the 
distributions  for  the  measured  temperature  fluctuation  and  the  measured  surface  height  are  very 
poor.  The  surface  height  is  the  classic  sinusoidal  cusp  shape;  whereas,  the  temperature  fluctuation 
appears  to  be  more  of  a  random  variable  distribution.  Clearly,  if  the  distributions  do  not  match,  it  is 
highly  unlikely  that  there  is  any  meaningful  relationship  between  the  surface  temperature 
fluctuation  and  the  surface  height. 

An  investigation  of  a  promising  technique  to  routinely  measure  ocean  surface  roughness 
characteristics  has  conclusively  demonstrated  that  shape  from  reflection  is  not  appropriate  for  fleet 
operations. 
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